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Abstract 



A theory for the many-electron multi-photon process is presented. It is shown 
that after single-electron excitation into some level in the continuum (ATI) an 
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, inelastic collision of the excited electron with the parent atomic particle can 



! result in an excitation of the ion. It may be the continuum state excitation 



giving the doubly charged ion or the discrete state which also greatly enhances 
the doubly charged ion production. The probability of these phenomena greatly 
exceeds that of the direct ionization of a single-charged ion. The single-electron 
ATI makes possible the two-electron process even in the moderate field. The 
example of two-electron excitations of He atom in a 780 nm laser field with 
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\ intensity above « lO^^W/cm^ is discussed 
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I. INTRODUCTION 



In this paper we consider many-electron processes in a low-frequency laser field. One of 
the most exciting manifestations among these phenomena is the creation of multiply charged 
ions in a laser field. Since its discovery (Suran and Zapesochnyi' 1975) it has been intensively 
studied experimentally for a variety of atoms and charge multiplicity of ions (L'Huillier et al 
1982,1983a,1983b, Lompre et al 1984, Luk et al 1983, Boyer et al 1984, Johann et al 1986, Chin 
et al 1985, Feldmann et al 1982a, Feldmann et al 1982b, Agostini and Petite 1984,1985a,1985b, 
Dexter et al 1985, Eichmann et al 1987, Camus et al 1989, Zhu et al 1989, DiMauro et al 
1988, Bondar' and Suran 1993, Walker et al 1994). Different experimental conditions such 
as a wavelength of radiation, intensity of laser field, duration of laser impulse were used. It 
was recognised that the ionization process strongly depends on both the atomic structure and 
the experimental conditions. Multiple ionization can proceed either as a sequential process or 
as a collective atomic response resulting in simultaneous removal of several atomic electrons. 
Noble gases are shown (L'HuiUier et al 1982, 1983a, 1983b, Lompre et al 1984, Walker et al 
1994) to exhibit simultaneous, direct removal of two electrons below the saturation level. After 
saturation the sequential ionization takes place (L'Huillier et al 1983a). The later statement 
agrees with the result of (Johann et al 1986) and with the theoretical estimations (Lambropoulos 
1985). For alkahne-earth atoms the two-electron mechanism has been recently shown (Bondar' 
and Suran 1993) to be predominant for doubly charged ion formation in the infrared range, 
while in the optical and ultraviolet ranges the sequential mechanism dominates (Feldmann et al 
1982a, Agostini and Petite 1985a, Dexter et al 1985, Eichmann et al 1987, Camus et al 1989, 
Zhu et al 1989, DiMauro et al 1988). 

The resonances can in some cases strongly manifest themselves in the formation of multiply 
charged ions. There is strong experimental evidence showing that creation of doubly charged 
alkaline-earth ions is strongly influenced by resonances in neutral atoms and single-charged ions 
(Feldmann et al 1982a, Feldmann et al 1982b, Agostini and Petite 1984, Dexter et al 1985, 
Eichmann et al 1987, Camus et al 1989, Zhu et al 1989). The calculation for 5- photon double- 
electron ionization of Ca atom are in agreement with the resonance mechanism (Tang and 
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Lambropoulos 1987). In contrast, the resonances seem to be much less important for ionization 
of noble gases where no resonance structure was reported (L'Huillier et al 1982,1983a,1983b, 
Lompre et al 1984, Walker et al 1994). 

There were several attempts to theoretically describe the collective mechanism of ionization 
in the absence of resonances. In (Boyer and Rhodes 1985, Szoke and Rhodes 1986) there was 
considered the possibility of energy transfer from the outer electron shell driven by the laser 
field to the inner-atomic shells. In (Kuchiev 1987) "the antenna" mechanism of ionization was 
suggested. According to this mechanism after single-electron ionization the ionized electron 
can absorb energy from the field and transfer it to the parent atomic particle in the inelastic 
collision. Recently the idea of rescattering of the first ionized electron on the ion has been 
discussed in (Corkum 1993). In (Ostrovsky and Telnov 1987a,1987b) there was developed an 
adiabatic theory of the multiphoton processes. In (Fittinghoff et al 1992) there was considered 
the possibility that the first ionized electron is taken away from the atom so quickly that the 
second electron is shaken off. Note, however, that this idea seems to contradict the adiabatic 
nature of the multiphoton ionization. 

In this paper a theoretical approach to describe the direct, collective many-electron process 
in the laser field is developed. The basic physical idea is the following. Suppose that first single 
electron ionization takes place. Suppose that then the ionized electron can undergo the inelastic 
collision with the parent atomic particle. Then one can expect the electron impact to result 
in the excitation of the ion into a discrete or continuum state. The reason for the excitation 
is the fact that the ionized electron strongly interacts with the laser field. It can absorb and 
accumulate high energy needed for ion excitation. The momentum of the ionized electron 
oscillates due to wigghng in the laser field. Therefore its "energy" is also time-dependent, it 
oscillates. The time-dependent "energy" of the ionised electron depends on two parameters. 
One of them is the time-dependent field momentum of the electron k = (F/u) sinuut, where 
F, (J are the strength and frequency of the laser field. The other is the constant translational 
momentum K. The total momentum is K -|- {¥/ uj) sin ut. Thus the kinetic energy of the 
photoelectron is Ej^(t) = (K + (F/u) sinc<jt)^/2. We will show that the ion excitation due to 
electron impact is possible if the ion excitation energy Eexc is equal or less than the maximal 
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kinetic energy of the photoelectron 




This is an important inequahty. It states that if the field is strong enough, F^/(2a;^) > E^xc, 
then we can forget about the translational momentum, considering K Ki 0. For weaker fields 
the translational momentum becomes vital. Here is a point where ATI, see (Agostini et al 
1979) and references in review (Freeman and Bucksbaum 1991), comes into play. ATI permits 
the photoelectron to occupy an excited level in the continuum. It is clear that the higher 
the level in the ATI spectrum is, the larger is the translational momentum. For sufficiently 
high level an inequality in Eq.(l) can be fulfilled permitting the mechanism to work. It is 
important that the necessary level in the ATI spectrum should not be as high as the ion 
excitation energy is, it can be well below the excitation energy. Really, the energy E'ati absorbed 
above the single-electron ionization threshold is equal to the averaged kinetic energy of the 
photoelectron, Eg^ti — -Br = + F^/(4a;^). The latter is always less than the maximal 

kinetic energy, £'k < max(£'K). The stronger the field is, the larger is the difference between 
£'k and max(£'K). As a result max(£^K) can exceed the ion excitation energy, while E^xi can 
be below it, E'ati < Eexc < max(£'K). All this means that absorption of low energy above 
the single-ionization threshold can make the considered mechanism to work resulting in high 
energy excitation of the ion. Section VII illustrates this statement using an example of double 
ionization of He atom. For the case considered absorption of only few (3-6) quanta above 
the single-ionization threshold is sufficient to allow for double ionization of He which needs 
very high energy. It is important that ATI takes place with high probability as the numerous 
experimental data presented in (Freeman and Bucksbaum 1991) and numerical calculations 
(Potvliege and Robin Shakeshaft 1989) demonstrate. This makes the considered mechanism to 
be very efficient. 

The idea of the process considered is illustrated by the Feynman diagrams in Fig.l. 
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The solid lines in Fig.l describe the behaviour of two electrons. First both of them are in 
the ground state of the atom. The lines marked a and b describe the corresponding atomic 
states. Then the electron a is ionised. The beginning of the ionization process is shown by the 
sloped dashed line representing the absorbed photon. The propagation of the ionized electron 
in the laser field, when it absorbs another photons, is shown by the double solid lines. Collision 
of the ionized electron with the parent atomic particle is shown by the vertical dashed line 
representing the Coulomb interaction between the electrons. The state marked c describes the 
ion excitation which can belong either to the discrete or to the continuum spectrum. In both 
cases it can strongly interact with the laser field. The diagram b describes the exchange process. 

The central question of the problem is whether this mechanism can work at all. It arises 
because in the process of ATI the ionized electron could go far outside the atomic particle 
and therefore inelastic collision with it could seem improbable. To prove that this mechanism 
works at all we have to fulfil an accurate calculation. It seems very complicated because we 
are to consider the perturbation theory over the electron interaction when the behaviour of the 
electrons is strongly influenced by the time-dependent laser fleld. Fortunately, the adiabatic 
nature of the problem results in a very important property of the two-electron amplitude: it 
may be presented as a product of the amplitude of single-electron ionization and the amplitude 
of electron ion impact in the presence of a laser fleld. Both latter amplitudes can be reliably 
calculated. As a result there appears the possibility of reliable ab inicio calculation of the 
two-electron process. 

Factorization of the two-electron amplitude is one of the most important results of the 
paper. The calculations proving it are fully presented because problems of the type considered 
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are not well-known in literature. We will be forced also to consider briefly two related problems: 
single-electron ionization and the scattering problem. 

Factorization of the amplitude makes it possible to find a reliable estimation for the prob- 
ability. It permits us to prove that the considered mechanism of two-electron ionization works 
well. There is a high probability for the first ionized electron to return to the parent atomic 
particle. The first ionized electron can be considered as a kind of antenna, it absorbs energy 
from the external field and transfers it to the ion. This physical idea was first considered in 
(Kuchiev 1987) but the corresponding calculations were not pubhshed. 

The electric field is assumed to be linearly polarized 



The frequency lo of the field is supposed to be low, uj « where / is the ionization potential. 
The Keldysh adiabaticity parameter 7 



is considered as arbitrary. The first ionised electron, which plays so crucial role in the consid- 
ered picture, will be described neglecting the static Coulomb field created by the parent atomic 
particle. A similar description was used in the Keldysh theory of single-charged ion formation 
(Keldysh 1965, Reiss 1980,1987) . One can consider it a reasonably good approximation when 
the non-resonant processes are important, for example for the two-electron ionization of no- 
ble gases. This approximation should be even more reliable for the problem of two-electron 
detachment from negative ions. 

In Section II the general formulae describing the single-electron ionization problem as well 
as the two-electron problem are presented. Section III contains the prove of factorization of 
the two-electron amplitude. Section IV presents the necessary consideration of the scattering 
problem. In Section V a qualitative physical discussion of the physical ideas governing the 
process of two-electron ionization is given. Section VI presents the estimation of the probability 
of the two-electron process for different regions of adiabatic parameter. ATI is shown to play 
the major role when 7 > 1. Section VII gives an illustrative example: the double ionization of 
He atom by a 780 nm laser field is briefly considered. 



¥{t) = Fcosa;^ . 



(2) 




(3) 
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II. MULTI-PHOTON PROCESSES 



In this section a formalism to describe the multi-photon processes when an atom is placed in 
a laser field is developed. The general expression for the amplitude Aji of an adiabatic process 
in the field is 

Af^ = ^J\^f{t)\Vpm^{t)) . (4) 

The wave function lipiit)) describes the initial state of the atom taking no account of the 
interaction of the atom with the laser field. This interaction is described by the operator Vplt) 

Vpit) = -rFcosut , (5) 

where r = is the dipole moment. We use the r-form for interaction with the laser field. 
The wave function (\E'j(t)| describes the final state of the reaction taking full account of the 
interaction of the atom with the laser field. It makes the amplitude (^) describe the multi- 
photon, nonlinear process. Eq.(^ looks similar to the usual expression for stationary processes 
except for the fact that in (^) there is the nontrivial integration over the period of time T = 
2n/LJ, where u is the frequency of the laser field. 

A. Single-electron ionization 

First let us consider shortly the well-known Keldysh problem of multi-photon single-electron 
ionization (Keldysh 1965). focusing our attention on those points which will be used in the 
following consideration of many-electron processes. The initial wave function for the single- 
electron process is taken as 

^i{t) = $„(r, t) = Mr) exp {-lEj) , (6) 

where </'a(r) is the wave function of the ground state of the atom, and Ea is the atomic binding 
energy. The wave function of the final state must include the interaction of the photo-electron 
with the external field. If one neglects the static ion field in the final state then this wave 
function is equal to the Volkov wave function (Volkov 1935) ^/(t) = $p(r,t), 
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<l>p(r,t) =exp|i {p + k^)r- Ep{T)dT + ^ |; 



(7) 



F 1 

kt = —smujt, Ep{T) = -{p + krf 
uj z 



which satisfies the Schrodinger equation 



^|-$p(r,t) = f-^A-rFcos^t)<l>p(r,t) . (9) 



at ' ^ V 2 

The additional phase pF/uj^ in Eq.(^ is chosen to provide a convenient property 

<l>p(r, t + T/2) = $-p(-r, t) exp{-zEpT/2}, (10) 

which is essential for consideration of the parity conservation law, see Eq.(|l^) below. Substi- 
tuting (|D,(0) in (^) one gets the known expression for the amplitude p) of single-electron 
ionization when I quanta are absorbed and the final-state momentum of the electron is p: 

^^^H^;P) = ^j\%{t)\VFma{t)) = ^£(p + k,|\/^(t)|0„)exp|l5M)} . (11) 

Here the usual definition |p) = expjzpr} is used and S{x) is 

/■^ 1 / F N 2 pF 

Six) = / -(pH sinx) dx — EaX . (12) 

JO 2 ^ a; ^ UJ 

The energy conservation law for the process of ionization reads 

E, + luj = Ep. (13) 

The parity conservation law for the ionization process manifests itself as a condition for the 
amplitude 

A(^)(/;-p) = (-l)'P,A(^)(/;p), (14) 

where Pa = ±1 is the parity of the atomic state a. It is easy to verify using Eq.(0) that the 
amplitude Eq.([TT|) satisfies this condition. 

The integrand in Eq.(|ll|) contains a large phase ~ 1/u; and therefore the steepest descent 
method is applicable. The saddle points, which must be taken in the upper semiplane of the 
complex plane x = ut, Im x > 0, in the region < Re x < 27r satisfy the equation 



S'{x) = l(p + -smxY -Ea = 



(15) 



which results in the following condition 



smx 



2 Nl/2 



(16) 



where p|| and p_L are the longitudinal and transverse components of the vector p in respect to 



the field F and k = ^21^^01. It follows from Eq.(|T5D that there are two saddle points xi,X2 in 
the region of interest which we will label in such a way that 



Recosxi > 0, Recosx2 < 0; Im xi,lm X2 > 0, < Re Xi,Re X2 < 27r . 



Calculating the integral in Eq.(|Tl]) by the saddle-point method one finds 



(17) 



(7 = 1.2 



(19) 



exp [-S{x,)\ {(p + kil^f (t)|0,)} 



where Al^\l,p),a = 1,2, is the contribution to the amplitude of the saddle point x = Xa^ 
satisfying Eqs.(|l5]), (|T7|) and is the part of the integration contour in the t plane which 
crosses the saddle point x^- (Note that evaluation of (p + kt\VF{t)\(j)a) needs accuracy because 
of the singular nature of this matrix element for complex t.) It is easy to verify that the 
amplitudes A^f\l;p) satisfy the condition 



A^2\l;-p) = i-iyPaAril;p) , 



I ID 



(20) 



which shows explicitly that the amplitude Eq.(0) obeys the parity conservation law Eq. ([T^) . 
The absolute values of A^i\l;p) and A'^\l;p) are equal. In contrast, their phases differ sub- 
stantially and this phase difference depends on the field F. In order to check it out consider 
small p, -C K^. For this case the phase difference is easily found to be 2 a/1 + •y'^Fp^^/u'^. The 
estimation for the smallest p\\ is p\\ ~ ^/uJ. Therefore the estimation for the phase difference is 
~ 2(F^/co'^ + K^)^/^/^/a; ^ 1. The large phase difference makes interference between A^\l; p) 



9 



and A2 {I; p) not important for many cases. Then calculating the probability one can suppose 
that 



Let us develop an approach to the two-electron processes in a way similar to the single- 
electron ionization. The main purpose of this section is to obtain the analytical expression 
describing the Feynman diagrams in Fig.l. The initial two-electron wave function may be 
chosen as a product of the single-electron wave functions 



This form of the wave function neglects the initial-state correlations, but they do not play a role. 
The most interesting for our purpose is the final-state interaction. To simplify the presentation 
we do not take into account explicitly the symmetry of the wave function with respect to 
permutation of electrons though it can be restored as discussed below. At the moment let us 
consider the two electrons as if they are distinguishable. The electron whose coordinates are 
ri will be considered as an atomic electron. It will also be referred to as "the first" electron. 
The electron with coordinates r2 will represent the degrees of freedom of the single-charged 
ion. It will be called "the second" electron. The wave function 0a(i'i) describes the ground 
state of the atom, Ea is the binding energy. The function (j)b{r2) is to be considered as the 
wave function of the ground state of the single-charged ion, Ef, being the binding energy of this 
state. Note that the removal of the first electron is an adiabatic process, see Section A. Hence 
the ion remains mainly in the ground state after single-electron ionization. We are interested 
in the interaction of the atom with the laser field. Therefore we consider the interaction of the 
first electron with the field. In particular we take into account the excitation of this electron 
into continuum states. In contrast, the direct interaction of the ion in the ground state with 
the laser field is strongly suppressed due to a high ionization potential of the ion. That is why 




(21) 



B. Many-electron processes 



ipi{ri, Fa, t) = 0a(ri)0fo(r2) exp {-i{Ea + Ei,)t)} . 



(22) 
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we will neglect direct interaction of the second electron with the laser field when it occupies the 
ground state. Certainly this interaction will be considered for the excited states. Thus the only 
possibility to be considered for the excitation of the second electron comes from its interaction 
with the first one. 

If we neglect the two-electron interaction then the final state wave function can be chosen 
as a product of the Volkov wave function for the photoelectron and the ion ground state wave 
function 

$/5(ri,r2,t) = $p(ri,t)0b(r2)exp{-zE;,t}. (23) 

Substituting Eqs.(p^),(p3l) into Eq.(|l) we reproduce Eq. ([TT|) for the single-electron process. 

Consider now the excitations of the single- charged ion into discrete or continuum states. The 
later case describes the doubly charged ion formation. The simplest wave function describing 
the excitation is given by the product of the wave function of the ionised electron and the wave 
function of the excited ion 

^^,(ri,r2,t) = <l>p(ri,t)<l>,(r2,t) . (24) 

Here we describe the wave function of the first electron in the continuum by the Volkov wave 
function (|^). The function $c(r2, t) describes the excited state c of the ion interacting with the 
laser field. It satisfies the time-dependent Schrodinger equation 

d 

i-^'^cir2,t) = {Hion{r2) - r2F cos tut) <l>c(r2, t) , (25) 

where -ffjon is the Hamiltonian of the single- charged ion. We will consider the discrete as well as 
the continuum state excitations. The wave function of the discrete state, if it is well separated 
from the other states, may be found in an adiabatic approximation 

$c(r2, t) = 0e(r2, t) exp j-z ^c(r)rfr| . (26) 

Here 0c(i"2,i) is the wave function satisfying the stationary Schrodinger equation 

Ec{t)(j)c{r2, t) = {Hion{r2) - rsF cos u;t)0c(r2, t) , (27) 
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with the time-dependent "energy" Ec{t). When the second order of perturbation theory over 
the field is apphcable then 

E,{t) =E,- ^ac{uj)F^ cos' ut , (28) 

where Ec is the energy of the stationary ion state c and a^uj) is the dynamical polarizability 
of this ion state. The excitation of the ion into the continuum will be described with the help 
of the Volkov wave functions (|^). It is convenient to consider the discrete state excitations 
and the excitations into the continuum on equal ground. Note that if we neglect the ion static 
field in Eq.(p7D then this equation is reduced to Eq.(|D describing the continuum states of the 
ion. Therefore we can consider the discrete states of the ion and its continuum states using the 
same set of Eqs.(p6|),(p7|) for the wave function and keeping in mind that ifion(i"2) ~ —(1/2)^2 
if a state in the continuum is considered. We will suppose that the phase in the definition of 
$c(r, t) is chosen in such a way that 

<l>e(r,t + T/2) = <l>p(,)(-r,t)exp{-z^,T/2}, (29) 

E,= ^ r E,it)dt, 
1 Jo 

where $ p(c) is the wave function of the state with the quantum numbers obtained by applying 
the inversion operator P to the quantum numbers of the state c. The definition Eq.(p9D agrees 
with Eq.(|l^). The possibility to satisfy Eq.(p9D follows from the Schrodinger Eq. (p5|) . 

The excitation of the ion takes place due to a mixing of the wave function ( P^ describing 
the excitation of the ion with the wave function (|2^) describing its ground state. In order 
to find this mixing consider the time-dependent Schrodinger equation for the final state wave 
function 

- ^-Q^mri, r2, t) = i^Hoir^, r^, t) + — j vl>* (n, r^, t) . (30) 

Here 

/fo(ri,r2,t) = -^^1 - riY cosujt + Hion{jC2) -r2Fcosu;t (31) 

is the Hamiltonian describing the considered two-electron system when the interaction between 
the electrons is neglected. This interaction is accounted for in Eq. (pOD by the term l/ri2. 
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Consider this interaction in the first order of perturbation theory. We are looking for the final 
state wave function which in the future, t oo, gives the excited ion described by the wave 
function Eq.fpll). Therefore the later wave function gives the zero approximation 

^;(ri,r2,t)^^f (ri,r2,t) = $p(ri,t)$,(r2,t) . (32) 

Then the correction 5\i/j(ri, r2, t) to the wave function caused by the interaction between the 
electrons satisfies the equation 

(^-^^ - ifo(ri,r2,t)^ 5M/}(ri,r2,t) = *(ri, r2, t) . (33) 

We are to find {S"^ f\ from Eq.(^) and use it in Eq.(0) for the amplitude. This procedure is 
simplified due to the fact that we neglect the interaction of the ion in the ground state with 
the laser field. This means that we neglect the transitions from the state described by the 
wave function 06 (r2). It makes the operator Vpit) in Eq.(|^) independent on the coordinate r2, 
Vp(t) = — FiFcostut. Thus the integration over r2 in (|^) is reduced to the projection of the 
final state wave function on 0b (r2) which is convenient to evaluate in advance. Let us denote 
it as ^}(ri,t) 

^}ir,,t) = 1 5^}(ri,r2,t)0fe(r2)exp{-^Ebt}rfV2 . (34) 
From Eq.(^) we find that it satisfies the equation 

(^-2^ + ^Ai + riFcosa;t^^}(ri,t) = K6(ri,t)$;(ri,t) , (35) 
where Kfe(i"i, t) is 

V,h{r^, t) = exp |z ^* ^,(r)rfr - E,t^ | j 0:(r2, t)^0fe(r2)rfV2 . (36) 

The solution of Eq.(^) reads 

mvr,t)= j^dt j %ir[,t)V,,{v[,t)Gir[,t;r,,t)dV, . (37) 

Here we introduce the time-dependent retarding Green function G{r[,t' ;ri,t) describing the 
electron propagation in the laser field. It satisfies the equations 
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-i— + -A + rFcostut G(r',t';r,t) = 6{r' -r)6{t -t), (38) 
cjij 2 J 

G(r',t';r,t) = 0, if t < t . 

Now from Eqs.(|37D,(^) we find the amplitude p, c) describing the two-electron process 

when n quanta are absorbed, one electron occupies the final state p in the continuum and the 
other occupies the discrete or continuum state c 

T oo 

A'^^'\n;p,c) = j;jdtj dt\<l>p{t')\V,b{t )Git ,t)VF{t)\<l>a{t)) . (39) 
t 

Integration over the variables ri,r[ is presented in (|39|) in the symbolic form, integration over 
the variable r2 is performed evaluating Vcb, see Eq. (|36|) . Note that one can restore the symmetry 
with respect to permutation of the electrons if along with the matrix element Vcb the exchanged 
matrix element is considered as shown in Fig.l b. 

The energy conservation law for the process considered reads 

Ea + Eb + nuj = Ep + Ec, (40) 

where n is the total number of absorbed quanta. Eq. (|40|) is similar to Eq.(|l^) for the single- 
electron process and can be justified in a similar manner, see Section III. 

Eq.(^) describes the ionization of the atom which is accompanied by the excitation of the 
ion into the discrete or continuum state. It gives the analytical expression for the diagram a in 
Fig.l. Let us formulate the main assumptions made in its derivation. 

1. The transitions from the ion ground state caused by the laser field are neglected. 

2. The wave functions in the continuum are described by the Volkov wave functions. 

3. Interaction between the electrons is considered in the first order of perturbation theory. 

III. FACTORIZATION OF THE AMPLITUDE 

Direct calculation of the amplitude ( ^9]) is complicated due to the adiabatic nature of the 
process under consideration which makes the integrand exhibit a strong variation with respect 
to the time variables t, t' . Fortunately the same adiabatic nature provides a very useful property 



14 



of the amplitude: the amphtude may be considered as a product of the single electron photo- 
ionization amplitude and the amplitude of electron-ion inelastic scattering. We will call this 
property factorization. For the first time this idea was considered in (Kuchiev 1987). 

In order to prove factorization let us first rewrite Eq. (^) in the more convenient form given 
below by Eqs. (pOf) , (|5TD , (|53D . With this purpose consider representation for the Green function 
through the Volkov wave functions 



G{r',t; r, t) = -te{t - t) J %{r , t>;(r, t) 



i2nY 



(41) 



This equality can easily be justified with the help of Eqs. (|38D ,([7|). Using it one can rewrite 
Eq.(ii 

A(2^)(n,p,c) = Jdt Jdt J y0^($p(t')lK.(t')|$,(O)(<f,(t)|Vi.(t)|$a(t)) . (42) 

i ^ ' 

It is important to distinguish in Eq. (^2|) the periodic functions of t' from the factors having 
the form exp{iEt'} which give the monotonic contribution to the phase. The function 0c(r2, t') 



defined in Eq.(^) is periodic. The last factor in the right-hand side of (|26D can be presented 



as 



exp |i J Ec{T)dT^ = exp {iEct} exp |z J {Ec{t) — Ec)dr 
Then the matrix element Vcfe(ri,t), see Eq. (p^) , can be presented as 



(43) 



V,b{v, i) = exp {i{E, - Eb)t}Ua{r, t ) , 



(44) 



where 



t/efe(ri, i) = exp <^ d / {E,{t) - E,)dT 



(>:(r2,t')— 0b(r2)dV2 
ri2 



is a periodic function of t! . Similar representation for the Volkov functions (0) reads 



'^'p(r, t) = Xp(r, t) exp {-iEpt} 



From 



Xp(r,t) = exp |z 
|),(||),(|3) we find 



[p + k,)r-|^* {Ep{ 



T - 



E^ rfr + ^ 



(45) 



(46) 
(47) 



15 



A(2^)(n;p,c) = i-t)^ Jdt Jdt J ^{x^(t')\U,,it)Mt')){xMVFma) x (48) 

t 

exp {t[{Ep + E,-Eb- E^)t + {E^ - Ea)t\} . 

Now it is useful to take into account the identity 

/ dt eyiv{iEt}f{t)=i 7f dtf{t)exp{t[{E-muj)t + mut]}- ^,(49) 

Jt T Jo E-muj + tO 



which is vahd for any periodic function f{t) = f{t + T). The identity can easily be checked out 
with the help of Fourier expansion f(t) = /m exp {— imcut}. Using Eg. (^91) we find from 
Eq.m 



A^'^\n;p,c)=J2A^'<^\nJ;p,c) , (50) 
I 

where 

^<...(„.,;p,e) = i- r,tdt7^fa"''l^'>"''l'^^"-''><^^,-'"lyf">l^°> X (51) 

K , J ^2 Jo J 27r3 Ep + Ec- Eb- E^- {n- l)uj + iO ^ ' 

exp {i[{Ep + Ec- Eb- Ea- {n- l)uj)t + {n - l)ujt]} . 

This formula looks similar to the usual representation for amplitudes of stationary processes. It 
is clear that the factor {xciif)\yF{t)\4'a) in Eg. (pl|) describes the matrix element of the transition 
of atomic electron into the continuum, while {Xvi^')Wc.b{^)\Xq{^')) is the matrix element of the 
inelastic scattering of the photoelectron on the ion. The non-stationary nature of the process 
manifests itself - there are nontrivial integrations over the time variables in Eg . (^l]) . The 
integration over t results in particular in the energy conservation law Ea + Eb + nuj = Ep + Ec, 
Eg . (^OD . To simplify notation in the following consideration the index of summation / in 
Egs.(0),(0) is chosen differently from Eg . (|49|) , I = n — m. The physical interpretation of the 
number / comes from integration over t which gives the /-th harmonic of the matrix element 
of single-electron ionization. Therefore I is to be considered as the number of guanta absorbed 
during the single-electron ionization. Similarly, the integration over t! shows that n — I is the 
number of guanta absorbed in the process of inelastic collision of the photoelectron with the 
ion. At this stage of consideration this interpretation is not of great importance, but later, 
when we prove factorization, it will become essential. 
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Our next step is to simplify Eg . (|5TD . With this purpose note the following. First, integration 
over t describes an adiabatic transition and therefore can be evaluated by the steepest descent 
method at the saddle points situated in the complex plane t. Second, the wave function Xqli"; t)^ 
see Eg ■ (^Tf) , depends on q in a very simple way. It is an exponent of a linear form of q, Xq(r, t) = 
exp{z[qR + a;]}, R and a being q-independent. This allows us to reduce integration over d^q to 
the residue at the pole arising due to a node of the denominator E^ = Ep + Ec — E}, — {n — l)uj. 
These facts make it convenient to deform the contour of integration over t into the upper 
semiplane of the complex plane t to bring it close to the saddle points 

j'^ ex^{ilujt}{xcSt)\VF{t)\(t)a)dt= I exp{ilujt}{x^{t)\VF{t)\(j)a)dt . (52) 

C2+C1+C0 

The integrand here is similar to one in the single-electron ionization problem, see Section II A. 
From its analysis we know that there are two saddle points of interest ut = xi,X2, see Eg . ([T7|) . 
We call Co-, a = 1,2 that part of the contour of integration which is situated in the vicinity 
of the point t = Xa/u. The integrand on the contour Ci is the function obtained by analytical 
continuation from those regions of the real axes t where cos cut > 0. Similarly, the integrand 
on the contour C2 is the analytical function obtained by continuation from those regions of the 
real axes where cos cut < 0. 

The contour Cq in Eg.(^) describes those parts of the contour of integration which are 
well separated from the saddle points and therefore give negligible contribution to the integral, 
/^^^ 0. There is no problem with the end-points t = 0, T as their contributions compensate for 
each other due to the periodic nature of the integrand. The saddle points ujt = Xi,X2 depend 
on q, but it is convenient to choose Ci, C2 as q-independent contours situated in the vicinity of 
the pole = Ep + Ec~ Eh — {n — l)uj. We can do it because then for any given q situated close 
to the pole we can deform Ci,C2 in such a way that these contours pass through the saddle 
points. As mentioned, the vicinity of the pole is the only important region of integration over 

q 

Substituting Eg.(0) in Eg. (|5T|) , using the coordinate representation given in Eg. (^T]) for 
the wave functions and integrating over d^q we find 
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1 

yl(2-)(n,/;p,c) = — J dt J dt J d'rdV x^ir' ,t)U,bir' ,t)i-rF cosut)^^) x (53) 

C1+C2 

i^ii? + lut + {n- l)ujt + \i r' -ktr- i -k^ 



27ri? 



exp < t 



ft 1 1 -pa 



where 



2 /cj 



1/2 



is the momentum of the photoelectron after absorption of / quanta. R 
function of all the variables of integration 

R = R(r, r'; t, t') = / k^rfr — r' + r = — (coscut — coscut') 



(54) 



in Eq.(^) is the 



r + r . 



(55) 



Eq. ( |53D is very convenient for evaluating parameters governing the process. First of all note that 
integration over the variables r, r' is localised in the vicinity of the atom. Really, integration 
over r describes the matrix elements responsible for the ionization process, while integration 
over r' describes the matrix elements of inelastic collision of the photoelectron with the ion. 
These r, r' should be compared with the amplitude of wiggling of the electron in the laser field 
F/uj"^. The latter becomes large even for quite moderate electric fields 

F 



> 1 



(56) 



Note that this inequality may be fulfilled for a big, 7 > 1, as well as for small , 7 < 1, adiabatic 
parameter. Therefore, for the fields satisfying (0) we can assume that 



3> r, r 



(57) 



This permits us to simplify R 



R ~ Rn — 



Rn ■ (v 



Rn 



Rq — Ro{t, t 



Ro = Ro(^, t ) = —(cos cut — cos cut ) 



(58) 
(59) 
(60) 



As a result the Green function in (^) becomes simpler 



i exp {zKiR} ^ exp {z (KiRo + • (r - r)) }, (61) 

= -K,^ . (62) 
/to 



Substituting Eq.(|6lD in (0) and using Eq.(p4D to rewrite the exponent we find 

T 

A(2-)(n,/;p,c) = j dtj dt {xp{t)\U,,{t)\Ki + 'kt)(}^i + \^t\VF{ma{t)) x (63) 

C1+C2 

r ' H "\ 

KiMt, t ) + + (Bp + E,~Ei- iA'?) (' - j' ikJdT I . 



Let us remember that integration over Ci is in the region where Re(cosci;t) > 0, while integration 
over C2 is in the region where Re(cosc(Ji(:) < 0. For inelastic collision the main contribution to 
the integral over t' comes from the vicinities of the points where 

I cos cut' I ^ , (64) 

see Eq. (^) in Section IV. (This equality has a clear physical meaning. If it is fulfilled then the 
wiggling energy of the photoelectron is maximal, F'^/{2uj'^), providing the best opportunity for 
the ion impact excitation.) Therefore we can assume Re(cosc<Ji(:) — cosujt' > when t G Ci, and 
Re(cosco't) — cos cut' < when t G C2. The function Ro(t,t') for complex t is to be considered 
as analytical continuation from the real axes. For real t the inequality Ro{t,t') > is fulfilled. 
As a result we find 

f? ( cos iut — cos iut\ if t G Ci , 
i?o(t,t) = — X (65) 

^ [ cos cut — coscijt, if t G C2 . 
Substituting Eq. (|65D in Eq. (|63D we come to the final expression for the amplitude 

A(2^)(n,/;p,c) = ^ E A^''\n,l;p,c), (66) 

a=l,2 

A'^^'\nJ;p,c)=A'^:\l-QaM'^'-''\n-l-Qa,uP,c)) • (67) 

Here a = 1,2 numerates contributions of the two contours Ci,C2 of integration over t. The 
quantities A'^^ are 

Ai^\i;Qa,i) = (<^>Q.,(t)|v>(t)|$.(t)) . (68) 



Ca 
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Comparing this expression with Eq.(|T9D we see that A^^^ are identical with the two terms in 
the amphtude of the single-electron ionization, see Eq.(|T8|). The momenta Qo-,« are defined in 
accordance with Eqs.(|0D,(|6|),(|6|) 

Q^,i = i-irKZ, a =1,2. (69) 

r 

The quantity Ri in Eq. (|67|) is equal to 

Evaluating this expression one can suppose that the non-exponential factor Ro(t,t') in Eq.(^) 
is Ro{t,t') ~ {— ly^^ {F/lu'^) cos Xa, cr = 1,2, see Eqs.(|6^), (|65D . The quantity cosXo- is to be 
found from Eq.(^) in which p\\ = {—lyKi, p± = 0, and the sign in the right-hand side is 
±1 = (—1)'^+^ in accordance with Eq.(|l^. As a result the non-exponential factor is the same 
for both saddle points: -Ro(''^ = Xi/uj,t') = Ro{t = X2/uj,t') = Ri] Re Ri > 0, Im Ri < 0. 

It follows from the definition of Ki ( [5^ ) that both amplitudes A^'' and A'2'^ in Eq. (|6^) 
describe the ionization process with absorption of / quanta. Note however, that in the considered 
two-electron problem they manifest themselves differently from the single-electron problem 
see Eq. ([T8|) . First, each of them depends on its own argument, which is the photoelectron 
momentum. It is Qo-,« for A^J^\ Qn = —Q2,i- Second, they appear in Eq.(|66D with different 
coefficients. 

The second factor in Eq.(^) is 

A(-2^)(n-/;Q.,,;p,c) = -^^ {%{t)\V,,{mci^,Mdt ■ (71) 

ZTT I Jo 

It can be recognized as the amplitude of electron impact on the single-charged ion in the 
presence of the laser field, see Eq.([78D in the next section. The initial and final momenta of the 
electron are Qo-,z and p. The collision results in the excitation of the ion into the state c and 
absorption of n — / quanta. The inequality ([56|) provides the following restriction on / 

1 



(72) 

J > 



_2cu3 V '2 

where [z]-^ is a minimal natural number greater or equal to z, z < [z]^. Eq.(|72D has a clear 
physical meaning: after absorption of I quanta the photoelectron must be in the continuum 
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state. Otherwise, if this inequahty is broken, the amphtude would be suppressed. There will 



appear the additional suppressing factor exp{ — y (/q — 1)ujF/uj'^} in the amplitude as it follows 
from Eq.(^). This factor is governed by the parameter 

:^ > 1 ■ (^3) 

The result of this section is the following representation for the amplitude of two-electron 
process 

A^'^\n- P, c) = E ^ E h P, c) . (74) 

l>l0 cr=l,2 

Here every term under summation over / describes the process when first / photons are absorbed 
during the single-electron photo-ionization process and then n — I quanta are absorbed during 
inelastic collision of photo-electron with the ion. The term with / = /q describes the process 
which starts from the near-threshold single-electron ionization. The other terms with I > Iq 
describe the processes when above-threshold single-electron ionization takes place. It will be 
shown in Section V that the role of the terms with / = Iq and those with I > Iq depends on the 
strength of the field: for a sufficiently strong field the term with / = Iq dominates, while for a 
weak field the ATI processes become extremely important. 

Every amplitude A^^^\n,l;p, c) in Eg . (fflD is presented in the factorization form: it is the 
product of the amplitudes of photo-ionization and inelastic collision, see Eq. (|67D . Due to large- 
ness of the phases of the amplitudes A^\n; Qo-,0 it is more convenient to take into consideration 
only one of them evaluating the other from the parity conservation law Eq.(^). Then Eq. (ff^) 
may be presented as 

Y: Af)(n,/;p,c)= (75) 

o-=l,2 

A^{\l; Qz, c) [A(^-2-)(n - I; Qr, p, c) + {-if PaA^'-^-\n - I; -Q,; p, c) = 
A^{\l; Qz, c) [A^--''\n - /; Q^; p, c) + {-1)^ PaPbA^'-''\n - I; Qi; -p. Pic)) . 

Here Qz = Qi,/ = —KiF/F, and Pb is the parity of the ion state b. The last line in Eq.([75D is 
evaluated taking account of Eq.(p]) which presents the parity conservation law for the inelastic 
scattering amplitude discussed in the next section. Eqs.([74D,(|75D are valid when the amplitude 
of the electron wiggling in the laser field is large, see Eqs.(^),(^). 
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IV. INELASTIC SCATTERING 



Consider the inelastic electron-ion impact in the laser filed restricting attention to the 
aspects of the problem which are most closely related to the many-electron excitations of the 
atom by the laser field. The general expression for the amplitude of inelastic collision in the 
first order of perturbation theory over electron-ion interaction is 



(76) 



where \E'j(t), \E'j(t) are the initial and final-state wave functions. Note that Eg . (|76D differs 
from Eq.(^) describing the general process. In Eq.([76D both and f{t) account for the 

interaction of the system with the laser field, while in Eq.(D only the wave function takes 
this interaction into account. The difference results from the fact that Eq . (|76D from the very 
beginning describes the perturbation theory over l/ri2. Note also that the factor — l/(27r) in 
Eq.(|76D is the usual one for a collision problem. 

The wave functions for the collision problem should be chosen as 



l^i(t)) = %{Ti,t)iPb{r2) exp{-iEbt} 
I*/W) = *p'(ri,t)^c(r2,t) . 



(77) 



Here p, p' are the momenta of the electron before and after collision and V'fe; are the wave 
functions of the ground and excited states of the ion. Substituting Eq.([77|) in ( [76D we find 

11/-^, 



A(-2^)(m,p;p',c) = --- / {^^'{t)\VM\%{t))dt 

Zn 1 Jo 



(78) 



1 



ri2 



p, (pbit)) exp { - Ssc{uJt)^ dt 



UJ 



where 



1 ^ ' F \' , , , 1 / F . V 

- p H smx + hJx/uj) p H smx 

2 \ to 2 \ UJ 



dx + ^P^P^. (79) 



UJ 



Eq. (|78D demonstrates that the previously considered in Eq. (|7TD quantity is the amplitude of 
inelastic collision indeed. The energy conservation law for the case under consideration reads 



2 '2 

— + Eb + muj = — + Ec , 
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where |m| is the number of quanta which are absorbed, if m > 0, or emitted, if m < 0, during 
the scattering. The parity conservation law for the inelastic scattering reads 

A^'^-^'\m, -p; -p', P(c)) = {-l)^P,A^--^'\m, p; p', c) , (81) 

where Ph is the parity of the state 6, and -P(c) is the state with the quantum numbers obtained 
by applying the inversion operator P to the quantum numbers of the state c. It is easy to verify 
using Eqs.(|10D,(|2|) that the amplitude Eq.(|78D satisfies Eq.(|81|). 



As is usual in adiabatic processes there appears the large phase ~ l/cj in the integrand in 
Eq. (|78D ■ Therefore the saddle points give the major contribution to t-integration. They satisfy 
the condition 

S'^^{ujt) = i (^p' + - sin ujt\ + Ec{t) - - (p + - sin cjt^ - Eb = . (82) 

It is fulfilled for those moments of t when the time-dependent "energy" of the system in the 
initial state coincides with the "energy" in the final state. If there is a real saddle point t then 
the process proceeds with a high probability. If it is in the complex plane, Im t > 0, then there 
is a strong suppression originated from the adiabatic nature of the process. This statement can 
be checked out for special cases when t-integration in Eq. ([78|) can be fulfilled analytically. One 
of them is the excitation into the discrete ion state c having small polarizability. Then one can 
neglect the time-dependence of the energy Ecii) and the corresponding wave function (pdt), 
Ec{t) ~ Ec, (f)c{t) ~ 0c and find from Eq.([78|) the simple representation 



A^^''^\m, p; p', c) = /(p; p', c)JUqF/cu')t^, (83) 

1 



/(p;p',c) = -^(p 



PAb) , (84) 



ri2 

where /(p; p', c) is the scattering amplitude in the absence of the laser field (considered outside 
the mass-shell), Jm{z) is a Bessel function and q = p — p' is the transferred momentum. 
Representation similar to Eq. (|83|) is well-known for the case of elastic scattering considered in 
(Bunkin and Fedorov 1966). A Bessel function is large when its argument is greater than its 
index. Thus the amplitude Eq.(p3D is not suppressed if |qF|/c<j^ > |^|- The latter condition is 
equivalent to one given in Eq.(^) for the considered case Edf) = Ec. 
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Consider now how much energy can be transferred to the ion from the electron and the field 
during collision. Suppose that the initial-state momentum p is given. The final-state moment 
p' will be considered as a free parameter and we will choose it to make the energy level Ec as 
high as possible. From Eqs.(|80|),(^) we find that to make the process probable the following 
inequality is to be fulfilled for some real moment of time t 

^{p + ktf >E,{t)-E,. (85) 

The quantity in the left-hand side of Eq. (^5|) may be considered as a time-dependent "energy" 
of the impact electron. To clarify the meaning of this equation consider first the discrete-state 
excitations of the ion. For the discrete state c we can use Eq. (|28|) to estimate its energy. Then 
from Eq.(^) we find the maximum energy gap between the ion ground state b and its excited 
state c which can be populated in the collision 



max{E, -E,) = l f + -) + ^pI • (86) 



2 

Evaluating Eq.(^6D it was assumed that aciuj)^'^ < 1. Note that in the left-hand side of 
Eq.(|86D there stands the non-shifted value of the energy Ec. The moments of time t for which 
the transition into the state c satisfying Eq.(|86|) takes place are those for which the wiggling 
energy is maximal 

|sino;t| = l. (87) 

Then we find that final-state momentum p' is to be 

F 



biil~— , p±^o 



Consider now the excitations of the ion into continuum state, i.e. creation of a doubly charged 
ion. For this case Ec{t) = (pi -|- kt)^/2, where pi is the momentum of the knock-out electron. 
It is interesting to find the maximum binding energy of the ion Ef, which would still allow the 
collision to result in double ionization. From Eq.(|85|) we find 

1 / F\ ^ 1 
max|Eb| = - i^pw + -) + ^pI , (89) 
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which agrees with Eq. (|86|) if we put in the latter one = 0. Conditions (|87D ,(|88[) are vahd for 
this case as well. Eg . (|87|) was used previously to evaluate the quantity Ri in Eq. ([70|) . 

Note that according to Eq. (|86D the maximum possible excitation energy of the ion is equal 
to the maximum of the "energy" of the impact electron. The latter is always higher, and 
can be much higher, than the averaged energy which stands in the energy conservation law: 
max(p + k()^/2 > Ep. This fact will be important in the following estimation of the probability 
of many-electron processes given in Section VI. 

V. QUALITATIVE DISCUSSION 

In the previous sections we considered the collective response of two atomic electrons to the 
laser field and found that the amplitude of the process may be presented as a product of the 
amplitude of single-electron ionization and the amplitude of inelastic scattering, see Eq.(^) or 
([75|). This means that at the beginning of the process single-electron ionization takes place. 
After that the ionized electron undergoes inelastic collision with the parent atomic particle. The 
general reason for factorization of the amplitude is quite clear. It follows from the existence of 
two different scales for electron localization. One is the radius of the atom. The other is the 
radius of the wiggling of the electron in the laser field which is supposed to greatly exceed the 
atomic radius. 

The simple physical meaning of the mathematical description presented above is the fol- 
lowing. The single-electron ionization may be considered as a quasiclassical movement of the 
atomic electron under the non-stationary, time-dependent barrier. When the tunneling is fin- 
ished the electron appears from under the barrier at some particular point separated from the 
atom. After that the electron propagates in a classically allowed region. It strongly interacts 
with the laser field exhibiting the wiggling. The propagation of the electron wave function in 
the classically allowed region may proceed in any direction because it starts from the point. In 
particular it propagates in the direction of the parent atomic particle. This makes possible the 
return of the electron to the single-charged ion. As a result the collision of this electron with 
the ion is very probable, the probability being proportional to the cross section and inversely 
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proportional to the square of the distance which separates the electron from the atom after 
finishing the under-barrier part of the trajectory. This distance is equal to the absolute value 
of Ri given in Eq.(|70|) 

\Ri\ = ^,{m + if+i']m-if+i']f": Pi = % 7 = 9^. (90) 

It depends on the amplitude of the electron wiggling F/ lJ^ as well as on 7 and the parameter /?/ 
which is the ratio of the photoelectron momentum to the field momentum. The photoelectron 
momentum depends on the number of absorbed quanta, which makes \Ri\ depend on it as well. 
This separation \Ri\ can exceed the amplitude FjiJ^ of the wiggling. For this case wiggling 
itself is not sufficient to prove that the electron can return to the parent atomic particle. That 
is why the discussion of all mathematical aspects of the problem given in the previous section 
is vital. 

There are two paths for the under-barrier electron escape from the atom. Roughly speaking, 
one of them is along the direction of the field F and the other is in the opposite direction. (The 
trajectory is described by complex time and coordinates. Therefore, it is more accurate to say 
that the real part of the trajectory is either along the field or in the opposite direction.) Every 
term A^^\a = 1,2 in the amplitude of single-electron ionization, see Eq.(|l^), describes the 
escape along one of these two paths. If we consider the pure single-electron ionization problem 
then after the tunneling the momentum of the electron should be the same for both routes of 
the escape. In contrast, for the problem of two-electron ionization the electron must return to 
the parent atomic particle and have the collision with it. The return route is in the direction 
opposite to the escape path, being different for the events described by A^i^ and for those 
described by A'^\ This explains why in Eq.(p7|) the momentum Qi,/, which is the argument of 
is opposite to the field direction while the momentum Q2,/, which is the argument of ^3'^^ 
is along this direction. 

VI. ESTIMATION OF THE PROBABILITY OF MANY-ELECTRON PROCESS 

There are different ways for the process of two-electron excitation with absorption of n 
quanta. During ionization of the first electron there may be absorbed some number / of quanta. 
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The remaining number of quanta n — / is absorbed in the coUision. The two-electron amphtude 
is the sum of the amphtudes of these independent events, see Eq. ([7^) . The interference of 
the terms with different / should be suppressed by the large and field-dependent phase of the 
amplitude Ai'\ see discussion in Section II A. Therefore as a first approximation we can neglect 
this interference. Then we find from Eq. ([74D the simple estimation for the probability W^'^^\c) 
of the two-electron process when one electron in the final state is in the continuum and the 
other is excited into the state c belonging to the discrete or continuum spectrum 

W^'<'\c) = J2W^''\n;p^,c), (91) 

n 

»•-'(„; P...) = E '-^('.0 = 0)^^^1^M^ . (92) 

l>lo 

Here the total probability W^'^'^\c) is presented as a sum of the probabilities with absorption of 
a given number n of quanta W^'^^\n;pn, c). The quantity dW^^\l, 6 = 0)/ d^l is the differential 
probability of the single-electron ionization with absorption of I quanta, and transition of the 
electron into the state with momentum Ki whose absolute value is given in Eq.(0), the condi- 
tion 6 = shows that the considered momentum direction is along the field, = KiF/F. The 
quantity cr(n — /; Ki;pn, c) in Eq.(p^) is the cross section of the electron-ion inelastic scattering 
in the presence of the laser field, Ki and Pn are the momenta of the electron in the initial and 
final states for the scattering problem. The scattering is accompanied by the absorption of 
n — I quanta and excitation of the ion into the state c. Evaluating Eq. (|92D we took into account 
Eq . (pID for the probability of the single-electron event and for the sake of simplicity neglected 
the interference of the two terms in the square brackets in Eq . (|75|) . The distance \Ri\ is given 
in Eq.(g). 

Consider now how high is the excitation energy provided by the considered mechanism. The 
general rules showing that the excitation is probable are given by Eqs.([86|),(^) in which the 
impact electron momentum for the considered case is = Ki, p± = 0. As a result we come 
to Eq.(l) 

1 / 

2[^^ + -) >^-c, (93) 
which it is convenient to present as 
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where the quantity 7exc niay be called the adiabatic parameter for the considered ion excitation. 
The energy E^^c is introduced to consider simultaneously the ion excitations into discrete and 
continuum spectra: i^exc = E^ — Ef, for discrete excitations and E^xc = 1-^61 for excitations into 
the continuum, compare Eqs.(^),(p9D. 

Eq. (^) demonstrates that the first ionized electron must possess "enough energy" to excite 
the ion. Due to wiggling in the laser field the "energy" of the electron is not conserved, it 
oscillates. Condition ( p3D states that at the most favourable moment of time, when the field 
momentum is big, kt = E/u, and parallel to the translational momentum Ki, this energy 
must be equal to the excitation energy. Note that the averaged energy of the photoelectron is 
always lower than the maximum one: {Ki + F/ujY/2 > Ej^^ = Kf /2 + F'^/{Auj'^). For further 
consideration it is convenient to distinguish the case of high wiggling energy from the case of 
low wiggling energy. Consider them in turn. 
1. High wiggling energy. 

Suppose that the energy of the electron wiggling is so high that it exceeds the ion excitation 
energy for the considered state c, F^/ {2lj'^) > E^ — Ef,, i.e. 

7exc < 1 . (95) 

Then it is sufficient to take into account only one first term / = /q in the sum over / in Eq. (p^) . 
It describes the process when first the near-threshold single-electron ionization takes place. The 
momentum of photoelectron is small Ki^ ^ 0. The electron wiggling provides enough energy 
to excite the ion during collision as is seen from the conditions (P^ , (^5]) . The probability of 
the process may be estimated as 

tyW(e)^^(;„.9 = 0) "<"-'°„^'^P-''> . (96) 

Eq . (|90D for the considered case reduces to 
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2. Low wiggling energy. 

Suppose that the wigghng energy is small compared with the excitation energy, F'^/2uj'^ < 
Ec — Eb, 

7exc > 1 • (98) 

Then the wiggling energy F'^/{2uj'^) is insufficient to provide the ion excitation. Therefore in 
this case we are to consider ATI during the single-electron ionization described by the terms 
I > Iq m the sum in Eg. (p^) . Let us call li the minimum number of quanta the absorption of 
which satisfies Eg. (p^) . From Egs.(p^),(^) we find 

F2 / , r 



> k 




J > 



7^ 



2^3 V 2 



(99) 



J > 



The ineguality (p^) shows that for a very weak field when 7exc >> 1, the above threshold energy 
must be egual to the ion excitation energy: ^ Ec — E^. In this case h can substantially 

exceed Iq- With increase of the field the necessary above threshold energy rapidly decreases. 
Really, it is seen from Eg.(^) that if 7exc ~ 1, then li ~ Iq. In this case the purpose of ATI 
is to supply the ionized electron with the momentum Ki which combined with the large field 
momentum will give the necessary high maximum "energy" for the photoelectron, see Eg. (p3|) . 
At the same time its averaged energy may be much lower than the ion excitation energy. The 
example of two-electron process in He considered in Section VII illustrates this possibility. We 
come to a conclusion that there are two reasons for increase of the probability of the two- 
electron process with increase of the field in the region Eg.(PBD. First, the total probability of 
single-electron ionization grows. Second, the stronger the field the lower becomes the necessary 
level in the ATI spectrum. 

The estimation for the probability of two-electron process is 

W'(-)(c) « ^Oi, 9 = O) "'" -''^^'^-'"'"'^' , (100) 

I ^1 I 

where Ri-^ is given in Eg. (pO]) . 

The total number of absorbed guanta n in the processes described by Egs.(p6|),( [10CI| ) must 
be large enough to provide the excitation of the level c as well as a sufficiently large final state 
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momentum p„ of the photoelectron: p„ ^ F/uo^ see Eg . (pSD . Therefore, the averaged energy of 
the scattered electron is -E'p,^ ~ "iF"^ / {Auj"^) . Estimation for the number of absorbed quanta n 
follows from the conservation law 



+ 



J > 



Ea — Ei,+ 



-I > 



2^3 l^^exc + 7+2 



J > 



(101) 



VII. DOUBLE-IONIZATION OF HE ATOM BY 780 NM LASER FIELD 

In order to illustrate the ideas discussed in the previous section consider the recent experi- 
mental results (Walker et al 1994) on precision measurement of double ionization of He atom 
by the 780 nm laser field in the intensity region from 0.15 to 5.0 PW/cm^. The data show 
the high ratio of He"*"^ to He"*" yield in this region of intensities. For the lowest intensity 0.15 
PW/cm^ the field momentum is F/uj = 1.2 atomic units. Therefore the maximal wiggling 
energy F^/{2uj^) = 19.3 eV is well below the lowest excitation energy of He+ which is 41 eV. 
Absorption of only 3 quanta above the single-electron ionization threshold changes the situation 
drastically. It gives the translational momentum K = 0.58. As a result the maximal energy of 
photoelectron becomes {K + F/oj)"^ /2 = 42.8 ev. The above consideration, see Eq.(p3D or ([l|), 
shows that this maximal energy of photoelectron can be used for ion excitation. Thus we see 
that absorption of only 3 above-threshold quanta during the single-electron ionization makes 
the excitation of He^ possible. The excitation might result in the double charged ion forma- 
tion due to field ionization from the excited state. Note that the above-threshold excitation 
energy itself in the considered example is sufficiently low. After absorption of 3 quanta it is 
only ~ 5 eV which is far below the excitation energy. This illustrates that the role of ATI is to 
supply the ionized electron with the momentum, rather than with above-threshold energy in 
ATI spectrum. Similarly, the absorption of 6 above-threshold quanta makes possible the direct 
excitation of the ion into the continuum, i.e. double-ionization. The probability of absorbing 
several quanta above the threshold is high for the considered region of intensities of the laser 
field, see for example the recent work (Schumacher et al 1994) where ATI for the first and sec- 
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ond harmonics of 1064 nm laser field with total intensity of 4 x 10^^ and 8 x 10^3 W/cm^ was 
examined. Thus the proposed mechanism qualitatively agrees with the experimental results. 

VIII. CONCLUSION 

We considered the mechanism of the many-electron process when first single-electron ioniza- 
tion takes place and then the ionized electron absorbs energy from the laser field and transfers 
it to the parent atomic particle. We proved by explicit calculations that this mechanism is 
very probable and provides the high energy necessary for ion excitation. The amplitude of the 
two-electron process is presented as a product of the amplitude of single-electron ionization 
and the amplitude of inelastic scattering. Both later quantities can be calculated ab inicio 
as discussed in the Sections II A and IV. Therefore the presented results provide a possibility 
for accurate calculations of many-electron processes in the laser field. The simple estimations 
for the probability of the process are also possible as discussed in Sections VI,VII. ATI proves 
to be very important for the mechanism considered. It supplies the ionised electron with an 
additional translational momentum. The latter gives sharp rise to the maximal energy of the 
photoelectron resulting in the possibility of excitation of the high-energy levels of the ion. 

The results obtained can be developed to cover related phenomena. One of them is ATI 
itself. We have seen that the scattering of the ionised electron on the parent atomic particle 
results in the transformation of its maximum energy into the excitation of the ion. Another 
possibility is to transform the maximum total energy into the energy of translational movement. 
This would strongly increase the probabihty of the population of high-energy ATI levels. For the 
first time this idea was discussed in (Kuchiev 1987). Another related problem is the high-energy 
gamma-quanta production by an atom in a laser field. 

In this paper we described the continuum states with the help of the Volkov functions. This 
description needs improvement, as the static Coulomb field in some cases is important, see for 
example (Reiss and Krainov 1994). One can expect this improvement to modify the amplitude 
of single-electron ionization and the scattering amplitude. It can also modify the parameter Ri 
describing the separation between the ion and first ionised electron. However, factorization as a 
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general property should remain valid. The reason for this comes from the fact that factorization 
is based on the parameter which is the large amplitude of classical wiggling of the electron in 
the laser field. For a sufficiently strong field this quantity can remain large enough even if one 
takes into account the ion Coulomb field. 

The other approximation made in the paper is the consideration of only the first order of 
perturbation theory over the electron interaction. This approximation can be insufficient if we 
are interested in the absorption of the very large number of quanta. The reason for this is that 
the electron impact transforms electron energy into the ion excitation, but there is always a 
possibility for the scattered electron to absorb more energy from the field and again transfer it 
to the ion during one more collision with it. 
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